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Background 

Ulcerative colitis (UC) and Crohn's disease (CD) are polygenic chronic inflammatory bowel 
diseases (IBD) of high prevalence that are associated with considerable morbidity. The 
hedgehog (HH) signalling pathway, which includes the transcription factor glioma-associated 
oncogene homolog 1 (GLIl), plays vital roles in gastrointestinal tract development, homeo- 
stasis, and malignancy. We identified a germline variation in GUI (within the IBD2 linkage 
region, 12q13) in patients with IBD. Since this IBD-associated variant encodes a GUI protein 
with reduced function and our expression studies demonstrated down-regulation of the HH 
response in IBD, we tested whether mice with reduced Glil activity demonstrate increased 
susceptibility to chemically induced colitis. 

Msthods and Findings 

Using a gene-wide haplotype-tagging approach, germline GUI variation was examined in 
three independent populations of IBD patients and healthy controls from Northern Europe 
(Scotland, England, and Sweden) totalling over 5,000 individuals. On log-likelihood analysis, 
GUI was associated with IBD, predominantly UC, in Scotland and England (p < 0.0001). A 
nonsynonymous SNP (rs2228226C^G), in exon 12 of GUI (Q1100E) was strongly implicated, 
with pooled odds ratio of 1.194 (confidence interval = 1.09-1. 3 1,p= 0.0002). GUI variants were 
tested in vitro for transcriptional activity in luciferase assays. Q1100E falls within a conserved 
motif near the C terminus of GUI; the variant GLI protein exhibited reduced transactivation 
function in vitro. In complementary expression studies, we noted the colonic HH response, 
including GUI, patched (PTCH), and hedgehog-interacting protein (HHIP), to be down- 
regulated in patients with UC. Finally, G//;^^'""^ mice were tested for susceptibility to dextran 
sodium sulphate (DSS)-induced colitis. Clinical response, histology, and expression of 
inflammatory cytokines and chemokines were recorded. G///+^''"^^ mice rapidly developed 
severe intestinal inflammation, with considerable morbidity and mortality compared with wild 
type. Local myeloid cells were shown to be direct targets of HH signals and cytokine expression 
studies revealed robust up-regulation of IL-12, IL-17, and IL-23 in this model. 

CORCiUSiOfiS 

HH signalling through GLIl is required for appropriate modulation of the intestinal response 
to acute inflammatory challenge. Reduced GLIl function predisposes to a heightened myeloid 
response to inflammatory stimuli, potentially leading to IBD. 



The Editors' Summary of this article follows the references. 
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Introduction 

Ulcerative colitis (UC; Mendelian Inheritance in Man [MIM] 
191390) and Crohn's disease (CD; MIM 266600) are chronic, 
relapsing, inflammatory bowel diseases (IBD) of high preva- 
lence (200-400 cases per 100,000 in Northern Europe and 
North America [1]) and are associated with considerable 
morbidity. Precise aetio-pathogenetic mechanisms are not 
understood but several lines of evidence implicate the central 
importance of a dysregulated host response to intestinal 
bacteria [2]. Epidemiological data, detailed molecular studies, 
and recent genome-wide association studies strongly suggest 
that UC and CD are related polygenic diseases that share some 
susceptibility loci (IL-23R, IL-12B, and NKX2.3). but differ at 
others: N0D2, ATG16L1, and IRGM are specific CD genes; the 
ECMl locus is associated with UC [3-10]. The IBD2 locus 
(Onhne Mendelian Inheritance in Man [OMIM] 601458) on 
Chromosome 12ql3 was first identified in a UK genome-wide 
linkage scan (peak LOD score 5.47 at D12S83) [11] involving 
both UC and CD patients. Later studies showed that IBD2 
contributes significantly to UC, notably extensive disease, but 
perhaps in a lesser way to CD susceptibility [12,13]. 

A strong candidate gene that maps to the IBD2 locus is 
GLIl, one of three related GLI transcription factors that 
transduce secreted hedgehog (HH) signals. HH signalling is 
key in gut development, homeostasis, and malignancy, but has 
not been carefully studied in IBD [14]. In developing 
intestine. Sonic (SHH) and Indian Hedgehog (IHH) provide 
a paracrine signal from epithelium to the mesenchymally 
expressed rei cpKn patched (PTCH) [14]. The binding of 111 1 
lo l'r(;il rc-lcases ihe membrane protein smoothened (SMO) 
from inhibition and allows HH signal transduction through 
the zinc finger transcription factors glioma-associated onco- 
gene homolog 1 (GLIl), GLI2, and GLI3 to direct tissue 
patterning and cell fate (Figure SI) (see [14] and citations 
therein). Chronic injury, inflammation, and repair are critical 
aspects of IBD, and thus it is pertinent that the HH pathway is 
centrally involved in these processes in several other tissues, 
including muscle [15], liver [16,17], and lung [18,19]. Indeed, 
HH signalling may play a central role in the inflammatory 
response since SHH is critical for T lymphocyte development 
[20], adult human CD4+ T cell activation [21,22], and myeloid 
ceU maturation in the spleen [23], Furthermore, SHH has 
recently been suggested to be a direct transcriptional target 
of NF-kB [24]. Dysregulation of components of the HH 
pathway has also been noted in inflammatory diseases of the 
gut, including Barrett oesophagus, chronic gastritis, and IBD 
[25]. Using microarray gene expression analyses of colono- 
scopic biopsies, we recently demonstrated that GLIl expres- 
sion is greater in the distal compared with the proximal 
colon, and is down-regulated in the intestinal mucosa in 
inflamed UC compared with noninflamed samples [26]. 

Taken together, these studies suggested a possible associ- 
ation between GLIl and IBD susceptibility, and we set out to 
directly test this possibility in both human disease and 
murine model systems. Our main initial objectives were to 
ascertain, from human genetic association study data, 
whether germhne GLIl variation was associated with IBD, 
and to describe the expression of HH signalling components 
in colonic inflammation in humans. These results led us to 
hypothesize that a reduced dosage of functional GLIl protein 
might play an important role in colonic inflammation. To test 



this hypothesis directly, we challenged GUI mice and 
their wild-type (WT) littermates with dextran sodium 
sulphate (DSS) to induce acute intestinal inflammation, and 
monitored clinical and histological parameters, cytokine 
profiles, and cellular targets of HH signalling in vivo. 

Methods 

Written, informed consent was obtained from all patients 
and controls. The study protocol was approved by the 
Lothian Research and Ethics Committee (LREC-2000/4/192 
and 2004/81103/22); Cambridge (LREC-01/418; MREC-03/5/ 
012); and Regional Ethics Committee, Karolinska Institutet. 

Participants and Samples 

Individuals for genotyping. In Scotland the population 
consisted of 817 IBD cases (474 UC; 335 CD; eight colonic IBD 
type unclassified PBDU]) and 1,374 healthy controls (HC) 
(Table 1). Diagnosis was made by clinical, radiological, and 
histopathologic al means, adhering to the criteria of Lennard- 
Jones [27]. All IBD patients attended the chnic at the Western 
General Hospital (Edinburgh, Scotland), a tertiary referral 
centre for IBD in southeast Scotland. The population was 
98.5% white, non-Jewish. The median age at diagnosis was 
31.0 y (IQR 23.3-46.0). Controls were recruited from across 
Scotland; median age was 50.0 (IQR 43.0-55.0), ethnicity 
>99% white, non-Jewish, and 48.7% male sex. 

In Cambridge, England 1,748 unrehiled Caiirasiaii IBD 
patients (Lennard-Jones criteria) of iiorllierii European 
origin altciidinglBD chnics in East ,\iiglia. UK were recruited 
(omprising 928 UC. 737 f:i). and 83 willi IBDU. The 
population was >99% white, non-Jewisli. I he 589 ethnically 
and geographically matched HCs were previously recruited in 
East Anglia for the European Prospective Investigation of 
Nutrition and Cancer (EPIC). Median age of controls was 60 y. 

In Sweden 288 UC and 205 CD patients (Lennard-Jones 
criteria) were recruited from various hospitals in Stockholm 
County. 281 Swedish HC were healthy volunteers from the 
Karolinska University Hospital staff (n = 170) and orthopae- 
dic day-case surgery patients (n = 111) with no previous 
medical conditions. Cases and controls were >99% white, 

Phenotyping. Phenotypic data were collected by patient 
questionnaire, interview, and case note review, and classified 
according to Montreal criteria (Table 1) [28]. 

Genotyping 

Scotland and Sweden. Genomic DNA was extracted from 
blood samples using a modified salting-out technique as 
previously described, and Nucleon kits. Genotypes were 
derived using the Taqman system for allelic discrimination; 
the assays were available from Applied Biosystems as Taqman 
SNP Genotyping Assays (7900HT sequence detection system; 
Apphed Biosystems), except for SNPs rsl0783819, rs3809114, 
rs507562, rs542278, rs730560, rsl 669296, and rs775322, which 
were genotyped on the lUumina platform. The accuracy of 
each Taqman assay was checked by repeat analysis in 5% of 
cases, with 100% concordance. Genotype distributions in all 
populations were consistent with Hardy-Weinberg equili- 
brium (HWE) (p > 0.01) for all SNPs. For the four tagging 
SNPs (tSNPs) that could not be derived by Taqman, genotypes 
were obtained by direct sequencing. 
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Table 1. Detailed Demographics and Phenotypic Data on Scottish, Cambridge, and Swedish IBD Population 



Cambridge Panel 



Swedish Panel 



Median age (y) at diagnosis/ 
recruitment (IQR) 



Ti (LI) 



■n (L2) 



lleocolon (L3) 
Upper Gl disease {LA 
Proctitis (El) 
Left-sided colitis (E2) 
Extensive colitis (E3) 
Inflammatory (B1) 
Stricturing (B2) 
Penetrating (B3) 



2,191 (1,374 HC; 474 UC; 2,337 (589 HC; 928 UC; 737 CD; 

Ipii CD; 8 IBDLI) 

HC 48.7%; UC 52.0%; CD 
39.6% 

HC 50.0 (43.0-55.0); UC 
34.1 (25.2 49.9); CD 27.8 

(20.8-41.1) liiiiiiiii;:; 

35.9% 
36.5% 



25.19 



HC 60.0 (53.0-69.0); UC 36.7 
126.84-50.35); CD 26.1 

i:;:i;2a3:-37.2) 
31.3% 
36.5% 
30.3% 



52.3% 
35.2% 
12.5% 



774 (281 HC; 288 UC; 205 CD; 
0 IBDU) 

HC 45.3%; UC 57.3%; CD 49.8% 



Gl, gastrointestinal; IBDU, colonic IBD type unclassified. 
doi:10.1371/joumal.pmed.0050239.t001 

Cambridge. DNA was extracted using Nucleon kits. 
Genotyping of CD cases and controls was performed using 
the Taqman biallelic discrimination system using an ABI 
7900HT analyser (Applied Biosystems). Genotyping of UC 
cases was performed using a 1,536-SNP Golden Gate bead 
array (Illumiiia). Concordance bclwccii pialfonns was as- 
sessed by genotyping 92 UC cases for SNI's rs222S22 1 and 
rs2228226 with concordaiue rales of 100% and 97.9% 
respectively and no evideiue of sysleinali< liias benveen 
platforms. Genotype distributions in case and control 
populations were consistent with HWE (p > 0.01) for all SNPs. 

Gene Expression by Microarray and Quantitative Reverse- 
Transcription PCR 

The cohort of patients used in the microarray studies 
consisted of 67 patients with UC, 53 with CD, and 31 HCs. For 
UC patients, eight were recruited at diagnosis (treatment 

naive), 18 had active disease, and 41 quiescent disease; ten 
were on systemic corticosteroid therapy, 11 on immunosup- 
pressants (azathioprine, mercaptopurine, or methotrexate), 
and 40 on 5-amino salicylic acid (5-ASA) therapy. Four of 53 
patients with CD were on corticosteroids, 13 of 53 on 
immunosuppressants, and 21 of 53 on a 5-ASA. Eight of 31 
HC had non-IBD inflammation: two scattered lymphoid 
aggregates with history of gastroenteritis, two microscopic 
colitis, one pseudomembranous colitis, one diverticulitis, one 
amoebiasis, and one eosinophilic infiltrate. 

The assessment of RNA quality and integrity was per- 
formed for the microarray experiments as follows: the 
amplified cRNA was purified using the RNeasy Mini Kit 
protocol (Qiagen) and 1 |il of amplified cRNA was quantified 
using the NanoDrop ND-1000 Spectrophotometer. The 
distribution of log intensities for each sample was plotted 
and outlier samples (i.e., greater than two standard deviations 
from the mean) were excluded from analysis. For additional 
demographics, phenotyping, and full methodology pertaining 
to the generation of the microarray dataset and quantitative 
reverse-transcription PCR (QPCR) see Noble et al [26]. 



Induction of DSS-lnduced Colitis and Histological Analysis 

Groups of two to four WT GliP^^ or WT littermate 
controls in mixed cages (bred at least three generations onto 
a C57BL/6 background) were administered 3% DSS in 
drinking water for 4-6 d. The amount of DSS consumed 
was not significantly different between WT and GliT^'^^ 
animals (unpublished data). Glil*"''^^ animals tended to weigh 
more than their WT littermates (mean. 28 g for GKi*^""^ 
animals and 24 g for \V 1 aiumals). llicrcforc. DSS-lrcated 
weight matched WT C57Bl/() aiumals {ii = \) were also tested 
along with Glil+^'"^ animals and WV Ijtlermates. No differ- 
ences were evident between WT littermates and weight 
matched C57B1/6. All animals were monitored daily for 
diarrhoea, bloody stool, and weight loss. Clinical scoring 
was as follows: 0, no symptoms; 1. diarrhoea; 2, bloody stool; 4, 
severe rectal bleeding and morbidity to the point of 
immobility/death. For histology, a segment of large intestine 
tissue of equal length and location for all animals was fixed 
overnight in 4% paraformaldehyde, dehydrated, infiltrated 
with paraffin, and sectioned at 5 |im. The remaining colonic 
tissue was taken for RNA analysis or frozen sectioning and 
Glil-LacZ expression analysis. Shdes were stained with 
haematoxylin/eosin and scored histologically by a gastro- 
intestinal pathologist (H.D.A.) blinded to the source of the 
tissue. Colitis scores were calculated using published methods 
[29]. 

Glil-LacZ Immunohistochemistry 

Colonic tissue was dissected out in cold PBS and fixed for 
30 min in 4% PFA at 4 °C. After a brief rinsing in PBS, the 
tissue was incubated in 30% sucrose/PBS at 4 °C overnight 
and then embedded in OCT and frozen on dry ice. 

Immunohistochemistry was performed on 8- to 10-)im 
frozen sections. Sections were dried, rinsed in PBS, and 
blocked for 1 h in 10% normal goat serum, 0.1% BSA, and 
0.3% Triton-X in TBS (50 mM Tris-HCl [pH 7.4], 150 mM 
NaCl). Antibodies used were: hamster anti-CD3 (Serotec, 
1:500), rat anti-CD19 (Serotec, 1:500), goat anti-CDllb 
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(Abeam 1:500), FITC-coiijugated anti-CDllc (Abeam 1:250), 
and rabbit anti-P-galaetosidase (a gift of J. Douglas Engel, 
Department of Cell and Developmental Biology, University of 
Michigan, 1:2,000). All antibody staining was performed over 
night at 4 °C, followed by incubation with fluorophore- 
labeled Alexafluor secondary antibodies from Molecular 
Probes, 1:1,000 for 1 h. Nuclei were counterstained with 
DAPl, and slides were mounted in Prolong Gold Antifade 
Reagant (Invitrogen). Images were taken using an Olympus 
BX51 at 200X for immunofluorescence and an Olympus 
FV500 confocal microscope at 1,000X for 0.30-^m optical 
sectioning. 

GLM Protein Mutagenesis, Luciferase Assay, and 
Immunofluorescence 

GLll EllOO was amplified from Image Clone number 
3531657, and cloned into pCMV-Tag2b, adding an N- 
terminal FLAG motif GLIl OllOO was obtained by inducing 
a point mutation using the OuikChange II Site-Directed 
Mutagenesis Kit (Siratagene) following the manufacturer's 
protocol. I'lasiuids encoding 8xGli-Luciferase, m8x Gli- 
Luciferase (mutated Gli sites), and Gli2AN were gifts of 
Andrzej Dlugosz (University of Michigan). 293T cells were 
plated in 12-well plates, transfected with 0.7 ^g/well tran- 
scription factor, 0.4 )ig/well reporter plasmid, and 2 ng/weU 
pRL-TK Renilla (Promega), and analyzed for luciferase 
expression 36 h after transfection using the Dual-Luciferase 
Reporter kit (Promega) following the manufacturer's proto- 
col. Firefly luciferase expression was normalized by well to 
Renilla, and fold changes were calculated by comparing to 
8xGli-Luciferase transfected alone. Statistical analysis was 
performed using the Student's f-test. For immunofluores- 
cence, 293T cells grown on collagen-coated coverslips were 
transfected with 1.5 )ig/well of transcription factor and 
analyzed after 36 h. FLAG-tagged GLIl isoforms were 
detected with mouse-anti-FLAG (Sigma) and anti-mouse 
AlexaFlour488 (Molecular Probes) both at 1:1,000. 

Cytokine Expression QPCR 

Colonic tissue was carefully cleaned of all fatty tissue and 
whole colonic mRNA was collected using Trizol, followed by 
RNA clean-up with DNase digestion using the RNeasy Mini 
Kit (Qiagen). cDNA was synthesized using the iScript cDNA 
synthesis kit (Biorad), and SybrGreen QPCR was performed 
on a Biorad iCycler. Expression levels were normalized to 
GAPDH, and statistical analysis was performed using the 
Student's Mest. 

Statistical Analysis 

Haplotype frequencies of the tSNPs were inferred using the 
expectation-maximization algorithm and these used to test 
whether haplotype frequencies were different in cases and 
controls as implemented in the EH and PM programmes. The 
test statistic 2 X (InLcase] + ln[Lcontrol] - ln[Lcase/Lcontrol]), 
which has a % distribution with « — 1 degrees of freedom 
(where n = number of possible haplotypes) was calculated and 
empirical j&-values obtained by permuting the data 10,000 
times. Haplotypes were examined and linkage disequilibrium 
between individual SNPs (r^) calculated using the Haploview 
programme v3.2 (http://www.hapmap.org). Individual SNP 
analysis was performed using % or Fisher's exact test, where 
appropriate, with two-tailed j&-values given and odds ratios 



(OR) presented widi 95% confidence intervals (CIs). The 
meta-analysis of SNP rs2228226 was performed using the 
Mantel-Haenszel method using a fixed effects model (R- 
software package). Details for calculation of false positive 
report probability (FPRP) are provided in Supporting 
Information (Text SI). Expression profiles were analysed 
using Mann-Whitney f7-test and Kruskal-Wallis test, assuming 
a nonparametric distribution of all datasets (GraphPad Prism 
4, GraphPad Software Inc.). 

Results 

Gene-Wide Variation in GUI Is Associated witli IBD and 
Attributable to a Nonsynonymous SNP (rs2228226) in tlie 
Scottisli Population 

Four multi-marker tSNPs (r^ > 0.8) were identified 
(rs38l7474, rs2228225, rs2228224, and rs2228226) to describe 
haplotypic variation of GLIl, detecting haplotypes of a 
frequency >1 %. We genotyped these four tSNPs in a Scottish 
IBD population consisting of 474 UC and 335 CD cases, and 
1,364 well-matched HCs (Table 1). We then used a model-free 
analysis [30] to test the association of GLIl and IBD 
susceptibility. In the Scottish population, we demonstrate a 
highly significant association in IBD {p < 0.0001) and UC (jb < 
0.0001), and an association with CD of borderline significance 
{p = 0.03). On analysis of individual estimated haplotype 
frequencies in Haploview, three common haplotypes were 
described (A to C; Table SI). We confirmed that this effect 
was confined to the GLIl gene by genotyping an additional 
seven haplotype-tagging SNPs, chosen from Phase II HapMap 
data, to tag neighbouring blocks, in 166 CD and 170 UC 
patienls. This confirnu-d llic |)rcsciu <• ol :i GUI S|);mning 
haplolyijc block lli:il did not cxlciid inio iu-iolibouriiio genes 
(LMIBE and AIUIGAP9) (Figure S2). 

The association on haplotype testing and log-likelihood 
analysis was largely attributable to a nonsynonymous SNP in 
exon 12 of GLIl (rs2228226C^G; tSNP4). rs2228226 was 
associated with IBD (allehc frequency OR= 1.23, CI 1.07-1.40, 

= 0.0026; homozygotes OR= 1.56, CI 1.15-2.11, p = 0.0047), 
CD (allehc frequency OR = 1.30, CI 1.08-1.55, p = 0.0053, 
homozygotes OR = 1.79, CI 1.21-2.65, p = 0.0048), and UC 
(allelic frequency OR = 1.19, CI 1.01-1.39, p = 0.04; 
homozygotes OR = 1.41, CI 0.98-2.03, p = 0.079) (Tables 2 
and S2). These data suggest an allele specific dose response 
with a greater OR for homozygotes than heterozygote 
patients. Mutation screening of the GLIl coding regions by 
direct sequencing failed to identify any novel SNPs. There 
was no association between seven additional GLIl variants 
from dbSNP and IBD (Table S2). Despite tight linkage 
disequilibrium (LD) across GLIl (Figure S2), there was notably 
less LD between rs2228226 and other tSNPs = 0.7) in both 
cases and controls (Figure S3), and this finding is likely to 
provide the explanation for the lack of association in 
neighbouring SNPs. 

Replication of GUI Association in Two Independent 
Northern European IBD Cohorts and Meta-Analysis 

We then sought to replicate these findings in other 
populations. In a large IBD panel from Cambridge, England 
(« = 928 UC, 737 CD, 83 IBDU, and 589 HC) association with 
GLIl was rephcated by log-likelihood analysis in IBD {p = 
0.009) and VC (p < 0.0001). rs2228226 was associated with 
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Table 2. Minor Allelic Frequencies for GUI Nonsynonymous SNP rs2228226 (tSNP4) in Scottish, English, and Swedish HC, IBD, CD, and 
UC 



Population HC IBD UC CD 

n Percent n Percent p-Value, OR (CI) n Percent p-Value, OR (CI) n Percent p-Value, OR (CI) 



Scotland 1,374 30.3 884 34.8 0X)026, TMWlMiM'M^ 0,042,1.19(1.01-1.39) 335 36.1 0.0053,1.30 (1.08-1.55) 

Cambridge 589 26.4 1,737 29.6 0.042,1.17 (1.00-1.36) 928 30.8 0.017,1.21 (1.03-1.42) 737 27.9 0.40,1.08 (0.90-1.28) 
Sweden 281 30.6 493 35.0 0.2i|lSt;;i:9§^j;i4S)i|ii5^^ 205 35.9 0.24,1.19 (0.90-1.58) 



Figure 1. Frequencies of estimated haplotypes in all three populations and thie full genotype data for the Scottish population are detailed in Supporting Information materials (Tables SI 
and S2). 

doi:1 0.1 371/joumal.pmed.0050239.t002 



IBD (OR 1.17, CI 1.00-1.36, p = 0.042) and UC (OR 1.21, CI 
1.03-1.42, p= 0.017), but not CD in this population (Table 2). 
As in Scotland, there was no association with tSNPsi_3. Minor 
allelic frequencies differed slightly from the Scottish pop- 
ulation; this difference is in keeping with that noted for a 
number of SNPs analysed for population stratification in the 
WTCCC study [3]. In the smaller Swedish cohort (n = 770), 
there was a consistent trend to association of rs2228226 with 
IBD (OR 1.14, CI 0.90-1.45) (Table 2) although this association 
was not statistically significant. However, the allelic frequen- 
( ics ill rases and controls were very similar in Scotland and 
Sweden iiulii ating that the lack of statistical significance in 
the latter cohort is due to insufficient statistical power. There 
were no consistent genotype-plienolype assoc iations across 
the populations, including slrali(i< alion for colonic LBl). 

A meta-analysis of Scottish, English, and Swedish data for 
rs2228226, using the Mantel-Haenszel method with a fixed 
effects model on IBD versus HCs (OR 1.194, CI 1.089-1.309, p 
= 0.0002) (Figure lA) confirmed the association of this variant 
in Northern Europe, achieving criteria for significance in a 
gene-centric study [31,32]. The meta-analysis was repeated 
separately for UC cases versus HCs with a similar effect size 
noted (OR 1.196, CI 1.077-1.327, = 0.0008) (Figure IB). 

Recognising the current problem with the publication of 
false positive findings in genetic association studies we 
estimated the probability that the association with disease 
risk found in the meta-analysis of GUI SNP rs2228226 
represents a true (rather than false positive) association by 
adopting the FPRP approach described by Wacholder [33], 
This gives an estimated probability that these findings 
represent a true finding of at least 92% (FPRP < 0.08). This 
method is designed to avoid over-interpretation of statisti- 
cally significant findings that are not likely to signify a true 
positive. In our study, this parameter gives clear support to 
our interpretation of these data. 

The GLM Variant Encoded by rs2228226 Is Functionally 
Deficient in Activating GLI-Responsive Transcription In 

Vitro 

rs2228226C^G is a missense mutation in exon 12 of GUI, 
encoding a change from glutamine to glutamic acid (QllOOE). 
The mutation falls within a well conserved motif at the C 
terminus of mammalian GTll proteins, near a recognized 
transactivation domain (Figure 2A) [34]. The QllOO residue is 
itself 100% conserved in all mammals examined (Figure 2B). 
In order to evaluate the functional consequences of the 



QllOOE mutation, we transfected either GLIl QllOO or the 
variant GLIl EllOO into 293T cells. No significant differences 
in level of expression or cellular localization were detected 
between these GLIl variants; both proteins were readily 
detectable in the nucleus of transfected cells (Figure 2C and 
2D). Blinded quantitation of cellular localization demonstra- 
ted a range of cytoplasmic and nuclear localization that was 
similar for both GLIl variants (unpublished data). We further 
evaluated the ability of each variant to activate the well- 
characterized GLI reporter SxGli-Luciferase [35]. We utilized 
Gli2AN, a very strong activator of Sx(;ii-Luciferase, as a 
positive control for GLIl transcriptional activity [36]. While 
both GLIl variants activated SxGli-Luciferase above baseline, 
WT GLIl QllOO was 2-fold more efficient as a transcriptional 
activator than the variant GLIl EllOO (Figure 2E). 

HH Pathway Activity Is Dysregulated in Colonic 
Inflammation 

Extended in silico analysis of microarray data demonstrates 
that mRNA transcripts of PTCH and hedgehog-interacting 
protein (HHIP), along with HH protein, are greater in the 
distal compared with the proximal colon in humans, mirror- 
ing the expression gradient of GLIl previously reported [26] 
(Figure 3A and 3B). GLIl, PTCH, and HHIP are pathway 
response elements whose expression levels predict pathway 
activity [14]. GLIl (;& = 0.0003). PTCH (p = 0.002), and HHIP (p 
= 0.0003) were lower in inflamed UC compared with HC from 
equivalent location (Figure 3C). IHH was lower in UC 
regardless of inflammation (p = 0.02). GLIl expression was 
lower in CD than HC (p = 0.004) irrespective of inflammatory 
status (Figure 3D), a noteworthy finding given that GUI 
variation was also associated with CD in Scodand. PTCH was 
lower in inflamed CD compared with noninflamed CD and HC 
(fi = 0.007) (Figure S5). GLIl and PTCH were both lower in 
non-IBD inflammation versus HC (Figures 3D and Figure S6). 
Treatment for IBD (5-amino salicylic acid, thiopurines, and 
corticosteroids) did not influence expression of GLIl (Figure 
S7). These data demonstrate overall down-regulation of HH 
pathway activity, including GLIl, PTCH, and HHIP, in areas of 
colonic inflammation. 

Qlj1+/iacz /^pjp^als Exhibit Mortality and Heightened 
Morbidity in Response to Intestinal Inflammation Induced 
by 3% DSS Treatment 

In vitro analysis of the GLIl llOOE variant demonstrated a 
50% deficiency in transactivation function compared to WT 
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Figure 1. Meta-Analysis of Nonsynonymous GUI SNP rs2228226 (tSNP4) 
in Scotland, Cambridge, and Sweden Using Mantel-Haenszel Method (n = 
5,352 Individuals) 

(A) IBD versus HCs. 

(B) UC versus HCs. There was no evidence of heterogeneity in the 

contribution of rs2228226 between the three cohorts (p = 0.825). 95% CIs 
for individual populations are represented by horizontal lines and 
population sizes by square boxes. The diamond represents the pooled 
OR (fixed effect model) with 95% CIs delineated by the diamond's width. 
Note the different range of the x-axis for (A) and (B). 
doi:10.1371/journal.pmed.0050239.g001 

GLIl and our genetic analysis demonstrated an allele-specific 
dosage response, suggesting that a moderate reduction in 
GLIl function may predispose to intestinal inflammatory 
disease. To specifically test this hypothesis, we treated 
Glil^^"'^ mice [37] and their WT httermates with 3% DSS to 
induce acute intestinal inflammation. 

Glil'^^'"^ animals were rapidly and severely affected by DSS 
treatment. After 6 d, four out of nine had died, and three of 



the survivors demonstrated severe morbidity, with significant 
rectal bleeding and almost complete immobility (Figure 4A). 
In contrast, no WT animals (n = 14) died, all were mobile and 
showed less morbidity on days 5 and 6 after treatment. 
Glil'^'"^^ animals developed bloody diarrhoea and significant 
weight loss by day 4, whereas WT animals did not develop 
clinical signs or measurable weight loss until day 6 (Figure 4B 
and 4C). The rapid onset of weight loss, clinical signs, and 
mortality in GHT^'^'^ animals was potentially indicative of 
extensive colonic inflammation, leading us to directly assess 
the histology of Glil^^ and WT colons after DSS admin- 

G//7"*"^'"^^ Animals Develop More Severe Colonic Pathology 
Than WT Littermates in Response to DSS Treatment 

Untreated WT and G/ii+"="=^ animals had intact colonic 
mucosa; no clinical or histological signs of inflammatory 
disease were present prior to DSS treatment (Figure 5A and 
5B). After 4 d of DSS treatment, WT colons exhibited 
evidence of inflammatory change but with few destructive 
lesions (Figure 5C), while extensive inflammatory infiltration 
and destructive colonic ulcers were prominent in Glil*"''^^ 
mice (Figure 5D). After 6 d, the number, size, and invasiveness 
of inflammatory lesions (Figure 5E-5G) were significantly 
greater in GliT^^ animals, and the overall colitis score 
(Figure 5H) [29] was significantly greater. Taken together, 
these results demonstrate that the loss of a single GUI allele 
leads to increased sensitivity to DSS treatment as reflected by 
severe intestinal inflammatory pathology and clinical signs. 

Intestinal Myeloid Cells Respond Directly to HH Signals 

We previously demonstrated that HH signalling is para- 
crine in murine intestine and colon [38]. Given that Glil*^^ 
animals contain a copy of bacterial P-galactosidase under the 
control of the GUI promoter, they represent an excellent 
model to confirm this observation and identify cell types that 
are Hh-responsive (i.e., LacZ positive) before and during 
inflammatory disease. These studies revealed that no auto- 
crine Hh signalling was detectable in untreated adult colon or 
following inflammatory stress; colonic epithelial cells do not 
express Glil or respond to Hh signals during colonic 
homeostasis or in areas of either mild or severe epithelial 
perturbation (Figure 6 and unpublished data). We did, 
liowever. detect a large population of Hh-responsive cells 
within the inflammatory infiltrate of DSS-induced ulcers after 
4 d of DSS treatment in Glil^^"'^ animals (Figure 6D-6F). 
Both lymphocytes [21.22] and myeloid cells [23] have recently 
been shown to respond to Hh signals outside of the 
gastrointestinal tract. We therefore investigated whether 
these populations were responding in the adult colon. We 
detected CD3-positive T lymphocytes, GDI lb-positive mye- 
loid cells, and CDllc-positive dendritic cells responding to 
Hh signalling in resting colon without DSS treatment (Figure 
6A-6C). After 4 d of DSS treatment, significant numbers of 
Hh-responsive myeloid and dendritic cells were seen in 
inflammatory ulcers (Figure 6E and 6F). We did not detect Hh 
response amongst the infiltrating T cells in DSS-induced 
ulcers (Figure 6D). Examination of CD19-postive B lympho- 
cytes in resting colon and after both 4 and 6 d of 3% DSS 
treatment revealed no clear Hh response in this population 
(unpublished data). Taken together, these data suggest that 
myeloid cells, including dendritic cells, are direct targets of 



PLoS Medicine | www.plosmedicine.org 



1766 



December 2008 | Volume 5 | Issue 12 | e239 



GLI1 and Inflammatory Bowel Disease 



Macaca molatta 
Mus musculus 
Rattus norvegicus 




C 




GUI nmi 




Figure 2. Q1100E Disrupts a Conserved Region of the GLI1 Protein and Reduces GLI1 Transcriptional Activity 

(A) Conservation of known functional domains in the Glil protein. Previously described Sufu binding, DNA binding, and transactivation domains 
[34,57,58] are shown schematically. Amino acid conservation of each domain is represented numerically and by shading of the bar below the domain. 
Red boxes indicate regions known to regulate GLI1 protein stability [49]. The conserved C-terminal domain that includes Q1 100E is adjacent to a known 
transactivation domain. 

(B) Alignment of the C terminus of mammalian Glil proteins. This region (AA 1080-1 106) is highly conserved in mammalian lineages. 
(C, D) GUI Q1100 and Ell 00 have similar cellular localization In 293T cells. 

(E) GUI E1100 is deficient in driving activation of the SxGli-Luciferase reporter compared to GLI1 Q1 100.GLI2AN is a strong activator of 8xGli-Luciferase 
and serves as a positive control for GUI activation. The mSxGli-Luciferase construct contains eight mutant Gli binding sites and serves as a negative 
control. Data are shown from six triplicate experiments done using two different plasmid preparations {n = 18). Bar heights indicate mean fold 
activation above baseline, and error bars indicate standard error of the mean. p-Values were calculated using the Student's t-test. 
doi:10.1371/journal.pmed.0050239.g002 



Hli signalling in the colon both during homeostasis and after 
inflammatory stress. This finding, in association with the 
findings above, implies that direct Hh signalling may play a 
role in innate immune function in the colon. 

G//7+'^'"'^ Animals Have Dramatically Increased IL-23p19 

and Pro-inflammatory Cytokine Expression after DSS 
Treatment 

Overall, unchallenged WT and Glil^^"'^ animals have a 
similar cytokine and chemokine profile as determined by 
QPCR examination of whole colonic tissue (Figure 7A). There 
are, however, several molecules (Cxcl2, Cxcl5, CxcllO, and 
Argl) that are modestly increased in G/ii+'*"'^ animals, and a 
small decrease in IT-6 expression; these changes may indicate 
a subtle difference in the inflammatory milieu of G/i2^'^'"^ 
animals at baseline. More strikingly, after 6 d of DSS Glil^'^^ 
animals demonstrate robustly up-regulated expression of 



several pro-inflammatory cytokines and chemokines com- 
pared to WT (Figure 7B). When baseline expression values are 
compared to values after DSS challenge for each genotype 
(Figure 7C), robust expression and up-regulation of ThI 
cytokines, including IFNy, is detectable in Glil'^^'^ animals 
(Figure 7B and 7C). We did not detect a significant difference 
in TGPP and IL-10 between Glil^^'^^ and WT animals, 
suggesting that down-regulation of anti-inflammatory cyto- 
kines was not the primary mechanism of increased inflam- 
mation in this model. 

The most highly expressed cytokine after DSS treatment in 
GUl^"'"^ animals was IL-23pl9, a molecule that drives 
differentiation of Th17 lymphocytes, key mediators of 
inflammation in several systems, including IBD [39,40]. IL- 
12 and IL-17, cytokines closely associated with IL-23, were 
also robusdy up-regulated in GliP''^'"^'^ animals. In fact, this 
entire pathway was highly up-regulated from the baseline in 
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Figure 3. Expression of HH Signalling Components in the Healthy Human Adult Colon (HC) and UC 

(A) PTCH, HHIP, and GUI mRNA levels increase along the length of the healthy adult colon, from ascending colon (AC) to descending colon (DC) and 
sigmoid colon (SC). 

(B) HH protein expression in terminally differentiated enterocytes at the luminal surface, is greater in the distal colon compared with the proximal 
(vertical and horizontal arrows represent these gradients). 

(C) Quantitative analysis of mRNA levels of IHH, PTCH, GLI1, and HHIP in UC compared with noninflamed HC (HC N-l). To account for the gradients 
identified along the length of the healthy colon (a-b), the data from SC only are shown. QPCR data is presented for IHH as this gene was not present on 

the Agilent microarray chip. 

(D) GUI expression in CD versus HC (N-l) and non-IBD inflammation (HC I) versus HC. Disease specimens are subcategorised into noninflamed (N-l) and 



PLoS Medicine | www.plosmedicine.org 



1768 



December 2008 | Volume 5 | Issue 12 | e239 



GLI1 and Inflammatory Bowel Disease 



inflamed (I) tissues. There was no change in levels of desert HH (DHH), PTCH2, GLI2, GU3, SUFU, or DISP1 In either UC or CD compared with HC, or in 
non-IBD inflammation (unpublished data). Analysis of SHH mRNA demonstrated a mild Increase In expression levels related to Inflammation that is 
consistent with the known expression of SHH In Inflammatory cells (Figure S4) [21]. Individual data points are plotted with horizontal lines representing 
the medians for each dataset. p-Values presented are derived from Kruskal-Wallls test, comparing levels In AC, DC, and SC, and from Mann-Whitney U- 
tests (UC versus HC (N-l)). 
doi:10.1371/journal.pmed.0050239.g003 



Glil^^"'^ animals compared to WT littermates (Figure 7B and 
7C). These findings are particularly interesting since the IL-23 
pathway has recently been strongly implicated in IBD 
pathogenesis both in humans [4] and mice [40]. Taken 
together, these cytokine data confirm that animals lacking 
one genomic copy of GUI develop significantly increased 
inflammatory disease compared to WT littermates in re- 
sponse to the same inflammatory challenge. 

Discussion 

The data presented here provide the first evidence, to our 
knowledge, that intact HH signalling is critical in the 
mammalian gut response to inflammatory challenge, and 
suggest that reduced GLIl function is implicated in IBD 
pathogenesis. We confirm that the HH signalling pathway is 
down-regulated in colonic inflammation in humans. We 
identify a specific GUI variant that is highly associated with 
UC/IBD, and demonstrate that the variant protein is func- 
tionally deficient as a transcriptional activator in vitro. 
Finally, we demonstrate that mice with a 50% reduction in 
GUI exhibit a heightened intestinal inflammatory response to 
DSS with significant up-regulation of the IL-23 pathway. Not 
only do these findings have clear implications for the 
understanding of IBD pathogenesis and suggest potential 
for therapeutic intervention, they are the first clear descrip- 
tion of a fuiu lional role (or III! signalling and C;i, 1 1 in bowel 
inflammation. 

The inherited variation in the GUI gene that we have 
detected is assotialed willi IBD and UC, in both Scotland and 
England, with findings lor rs2228226 confirmed by meta- 
analysis of over 5,000 individuals, with OR of 1.19. Evidence 
for an effect in CD is seen in the present study, but the 
predominant effect is clearly related to UC, consistent with 



the reported IBD2 linkage studies [11,13]. The magnitude of 
this association is entirely in line with the effect size noted in 
a number of recent studies of complex disease genetics, 
including CD [41], colo-rectal cancer [42], and coeliac disease 
[43]. The level of significance attained satisfies suggested 
criteria of p < 10"^ to 10"® for gene-centric studies [31,32]. 
Population heterogeneity has been previously described at 
the IBD2 linkage region [11,13,44,45]. The three Northern 
European populations studied here have previously demon- 
strated similar contribution of other IBD susceptibility genes/ 
loci, including NOD2 and IBD5 [46]. 

Whilst our resequencing efforts identify rs2228226 as the 
only coding variant associated with IBD, further deep 
resequencing will be important in the future as the haplotype 
analysis and log-likelihood analyses raise the possibility that 
other germ-line variants may also contribute to IBD risk. 
These need be explored formally — specifically the role of 
intronic variants, long-range promoter effects, and/or copy 
number variation. In this context, several complex disease 
genes, including NOD2 [8,9], have multiple independent 
mutations conferring disease risk, some disease genes have no 
causative mutations within coding sequences (e.g., IRGM in 
CD [10,47]), and synonymous SNPs may be associated with 
functional effects [48]. 

rs2228226C^G encodes a change from glutamine to 
glutamic acid (QllOOE). Our in vitro data demonstrate that 
GLIl llOOE is a subfunctional transcriptional activator 
compared to WT GLIl, though it is appropriately synthesized 
and localized. The QllOOE mutation causes a significant 

known tranLctivation region of GLIl; this change could 
directly modify transactivation activity, disrupt the structure 
of the transactivation domain, or affect protein stabilization 
[49], decreasing activity. Our murine data confirm that 
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Figure 4. GliV^'"'^ Animals Show Mortality, Severe Clinical Symptoms, and Profound Weight Loss after DSS Treatment 



(A) Kaplan-Meier survival curve. WT animals a 
in response to 3% DSS treatment for 6 d. 

(B) Glil^^'°^ animals display markedly more s 
severe bleeding/death. Each dot represents 

(C) Glil^^'"^^ animals (n = 9) lose weight moi 
doi:10.1371/journal.pmed.0050239.g004 



? 100% viable over the 6-d treatment period (n = 14). Nearly 50% of G//7+'''''^^ animals (four out of nine) die 

vere symptoms than WT animals after 4 or 6 d of 3% DSS treatment. 1, diarrhoea; 2, bloody diarrhoea; 4, 
n individual animal and the solid line shows the mean observation in each cohort. 
' rapidly than their WT llttemnates (n = 10). *, p < 0.05 by Student's f-test. 
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Figure 5. Glil^^'"'^^ Animals Demonstrate More Severe Intestinal Inflammation than WT Littermates in Response to DSS Treatment 
(A, B) WT and G//7+'''°^^ animals demonstrate normal thickness and structure of colonic mucosa in the absence of DSS. 

(C) WT animals (n = 4) exhibit mild colonic inflammation but do not develop substantial epithelial or ulcerative inflammatory pathology within 4 d of 
DSS treatment. 

(D) G//7+'''°^^ animals (n = 4) develop significant inflammatory infiltration, epithelial damage, and ulceration within 4 d of DSS treatment. 

(E) WT animals (n = 14) demonstrate ulceration following DSS treatment. 

(F) Glil^^'°^ animals develop profound intestinal inflammation in response to 3% DSS treatment, with severe epithelial damage in long stretches of their 
colonic mucosa (n = 9). 

(G) Blinded histological scoring of colonic damage after 6 d of DSS treatment. Standard lengths of tissue from the mid colon and distal descending 
colon were scored in each animal. G//!^'''''^^ animals (n = 6) have more overall inflammatory foci and more long foci (10+ crypt units affected) than WT 
animals (n = 6). Each dot represents the number of observed foci in an individual animal; the solid line shows the mean observation in each cohort. Red 
dots indicate the animals that were analyzed for cytokine expression. 

(H) Comprehensive colitis scoring [29] of WT and Glil'^'"'^ animals after 4 and 6 d of DSS treatment. p-Values for (G) and (H) calculated by the Student's 
doi:10.1371/journal.pmed.0050239.g005 
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Figure 6. Intestinal Myeloid Cells Respond Directly to Hii Signalling during Homeostasis and Inflammation 
(A-C) No DSS treatment. 
(D-F) 4 d of 3% DSS. 

(A, D) Few CD3+ T lymphocytes respond to Hh signalling during colonic homeostasis. Infiltrating T cells in DSS-induced ulcers do not express Glil or 
respond to Hh signalling. 

(B, E) Many but not all CD11b+ myeloid cells respond to Hh signalling and express LacZ. After 4 d of DSS, ulcers in Glif^"'^ animals are filled with 
myeloid cells, many of which respond to Hh signals. 

(C, F) CD1 1c+ dendritic cells respond directly to Hh signals during homeostasis and inflammation. Nuclear DAPI staining was performed but has been 

omitted to allow clear visualization of co-expression of nuclear LacZ and cell surface markers. 

doi:10.1371/journal.pmed.0050239.g006 



reduction in GLIl activity has significant consequences for 
the inflammatory response. 

Intestinal HH signalling in mice and humans is exclusively 
paracrine: HH ligands in the epithelium signal to the HH 
response network (PTCH, GLIl, HHIP) m the mesenchyme 
[38,50,51]. This has been demonstrated during development 
[38], homeostasis (Koltemd A, Crosse AS, Zacharias WJ, 



Walton KD, and Gumucio DL, personal correspondence), 
malignancy [51], and now, in this report, in inflammation. 
The human expression data presented here also support the 
relevance of the HH pathway to colonic inflammation in 
humans. GLIl expression, a reliable indicator of pathway 
activity, IS down-regulated in all forms of colonic inflamma- 
tion examined (UC, CD, and non-IBD inflammation). Whilst 
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Figure 7. Cytokine Analysis of GliV '"'^ and WT Mice after DSS Treatments Demonstrates Robust Pro-inflammatory Cytokine Activation 
(A, B) Cytokine and chemokine expression normalized to GAPDH is plotted on the y-axis for GM^^'"^ mice, and on the x-axis for WT animals. Gene 
expression levels that are statistically different are shown with stars. The dotted diagonal trend line indicates identical expression levels between WT 
and GhT^'"''^ mice. (A) Cytokine and chemokine expression in WT and GU'\'^''°'^ animals (n — 3) without DSS treatment. Cytokine and chemokine 

expression is low and very similar between WT and Glil^^'"^^ animals. 

(B) Cytokine expression in WT and G//7+'''°^^ animals (n = 4) after 6 d of 3% DSS treatment. Many pro-inflammatory cytokines and chemokines are 
significantly more highly expressed in G//7+'''°^^ animals compared to WT. Note the difference in axes in panels (A) and (B). 

(C) Fold change of cytokine expression from baseline to inflamed for WT and Glil^'"'^^ animals. We detect dramatic up-regulation of Th17 and Thi 
pathway cytokines and chemokines in Glil^^'"^^ animals. *, p < 0.05; **, p < 0.01 by the Student's f-test. 
doi:10.1371/journal.pmed.0050239.g007 
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embryonic HH pathway activity is reportedly recapitulated in 
diverse organs/tissues in response to acute injury and 
inflammatory challenge [15-19], GLIl expression is notably 
decreased in psoriasis, a disease of the epithelial border that 
shares some pathogenic features with IBD [52]. 

The localisation of HH protein demonstrated here is 
similar to that shown for IHH by van den Brink and 
colleagues by in situ hybridisation [53]. However, along with 
other commentators, we would urge caution in the use of the 
commercially available HH pathway antibodies presendy 
available [54]. We note further that the expression analysis 
of disease samples in humans is inherendy limited in that it 
only depicts a complex process at a late time point in 
established disease. Nonetheless, in combination with studies 
utilizing models (such as the GliP*^^ mouse) to directly 
correlate HH response with inflammatory challenge, these 
studies provide useful information to help elucidate the 
aetiology of human inflammatory disease. 

Glil^^"'^ animals, which have only 50% of the WT level of 
Glil, develop severe and rapid inflammation in response to 
DSS, despite having normal histological structure and a 
largely WT cytokine and chemokine profile prior to DSS 
treatment. These data, to our knowledge, provide the first 
description of a phenotype for reduced Glil function [37], and 
demonstrate the key role that a full complement of Glil plays 
in protection from inflammatory disease. Our in vitro data 
demonstrate that GLIl EllOO is capable of activating some 
GLIl response, suggesting that under homeostatic conditions, 
GLIl EllOO could perform adequately. Similar to the situation 
in GliT^^'^^ animals, however, under conditions of inflamma- 
tory stress, GLIl El 100 can only function at 50% of the level of 
WT GLIl. Whether the predisposition to inflammation in 
these systems is a direct result of lowered HH signal 
transduction within HH-responsive inflammatory cells or 
reflects the effect of lower HH signals on other stromal target 
cells [15] that, in turn, release signals that impact the integrity 
of the epithelial layer or the inflammatory phenotype of 
resident immune cells is not yet clear, and will be an 
important area of future study. Furthermore, we believe that 
it will be important to determine if other primary (e.g., germ- 
line variation in other HH pathway genes and the HH 
interactome, epigenetic changes, and gene-environmenlal 
interactions) or secondary phenomena (e.g., due to alterations 
in commensal flora) in the human/murine intestine alter the 
transduction of HH signals in response to acute and chronic 
inflammatory challenge. 

We have additionally shown that the HH pathway may 
directly modify the innate immune response through signal- 
ling to myeloid target cells. This finding is in actcird willi 
recent data demonstrating a crucial role for Hh signalling in 
myeloid cell maturation in the spleen [23]. Interestingly, 
myeloid cell populations are known to modify the intestinal 
inflammatory milieu through a significant impact on the IL- 
23/IL-17 pathway [55], a pathway that we find to be clearly 
and distincdy dysregulated in G/ii^'^"'^^ inflammation. Taken 
together, these data raise the exciting possibility that the 
normal role of HH signalling is to promote a tolerogenic 
phenotype in mucosal myeloid cells. In the face of reduced 
HH signal transduction in these cells, minor inflammatory 
stimuli may trigger robust inflammatory responses. Future 
studies directly investigating this hypothesis by specifically 
testing the effects of Hh signalling on intestinal myeloid cells 



will be crucial to the evolution of a more complete under- 
standing the role of Hh in intestinal inflammation. 

Limitations of the Study 

A number of issues meriting further work emerge from the 
pn^senl study. Firsdy, deep resequencing of the entire GLIl 
gene is rt^quired to identify additional variants that contrib- 
iilc lo Itu' overall disease association noted for the gene. 
Secondly, the in vivo examination of reduced Glil function is 
presently limited to heterozygous mice. Despite the robust 
phenotype observed in these mice following acute inflamma- 
tory challenge, these findings need to be extended to 
homozygous null mutant animals (GUr'^). Thirdly, additional 
animal models of colitis, particularly of chronic disease, need 
to be examined; for example, trinitrobenzene sulfonic acid 
(TNBS) induced colitis, or adoptive transfer of 
CD4+CD45RB'''S'' T cells into recombinase-activating gene 
(rag)-\ deficient (ragr'^) mice. Finally, we have not yet directly 
examined the phenotype of reduced GLIl function in antigen 
presenting cells in vitro. 

Conclusions 

A key paradigm emerging in complex disease genetics is 
that gene discovery may lead to fundamental shifts in the 
direction of basic and translational research efforts. Recently, 
it is apparent that mechanistic insights may have greater 
significance than the index gene discovery itself, especially 
when one considers the small contribution each locus/variant 
makes towards disease susceptibility. IBD has provided 
several recent examples. Notably, the discoveries of N0D2 
and IL-23R have launched academic and pharmaceutical 
research efforts into innate immunity and Thl7 signalling, 
whereas ATG16L1 and IRGM have focused much attention on 
the role of autophagy in CD pathogenesis. In the present 
study, our genetic studies led in turn to detailed observations 
of expression in established human disease as well as a 
mechanistic study in the mouse. We believe that these 
combined findings have the potential to open entirely novel 
lines of enquiry that may have great consequences for the 
understanding of IBD pathogenesis and potentially for 
inflammatory pathways in other HH-expressing organs (e.g., 
skin. lung. bone). 

In conclusion, we demonstrate here llic dvsrcgiihilioii of a 
developmental signalling pathwav in IBD lor llic lirsi lime to 
our knowledge. Furthermore, we show that these effects are in 
part genetically determined, with evidence implicating GLIl 
as an IBD2 gene, and identification of a specific variant with 
rt^diicc^d Iraiiscriptional activity. Itic Itmctional relevance of 
GUI is demonstrated by the severe intestinal inflammation 
that develops in the face of a 50% reduction in Glil 
concentration in an established mouse model of colitis. Taken 
together, these data strongly argue that tolerance to inflam- 
matory stimuli requires a fully functional HH signal trans- 
duction network. With HH signalling modulators under active 
investigation in a number of clinical trials (http://www.curis. 
com and http://www.clinicaltrials.gov), these data provide 
promise for the development of novel therapeutic strategies. 
Since the use of HH inhibitors has recently been advocated for 
the treatment of certain tumours that secrete HH, our finding 
that reduced HH signalling has inflammatory consequences 
suggests that inflammatory status should be carefully followed 
during upcoming Phase II trials [56]. 
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Supporting Information 



Figure SI. Schematic of Activation Events in tiie HH Signalling 

There are three mammalian homologues of Drosophila gene Hh. of 
which SHH and IHH have important functions in the gut and the 
immune system. Intestinal HH expression is limited to epithelial cells, 
while HH-response genes are expressed only in intestinal mesen- 
chyme. In epithelial cells, soluble HH is processed in the Golgi and 
released via Dispatched (DISPl). In the mesenchyme, upon HH 
binding to the transmembrane receptor PTCH, PTCH can no longer 
repress SMO, and SMO prevents cleavage of transcription factors 
GLI2 and G03 from activator to repressor form; GU activator forms 
then predominate, and transcription is activated in a raft of 
downstream targets. Direct early targets of HH signals include 
GLIl, a strong transcriptional activator which acts an amplifier of HH 
signalling PTCH, and HHIP, a feedback inhibitor; these markers can 
be used as indicators of pathway activity [14,54]. 
Found at doi:10.137l/journal.pmed.0050239.sg001 (66 KB PPT). 

Figure S2. The Haplotype Structure in GLIl and Surrounding 
Haplotype Blocks Are Presented in Haploview from Genotyping in 
UC (A) and CD (B) 

There is a haplotype block spanning the GLIl gene that does not 
extend into neighbouring genes. LD is described here by D' values. 
Found at doi:10.137l/journal.pmed.0050239.sg002 (1 MB TIF). 

Figure S3. The Haplotype Structure of the Four GUI tSNPs from the 
Scottish HCs and UC Population and from HapM;!]) CKL Is 
Demonstrated with Values 

The LD between rs2228226 and neighbouring SNPs is overestimated 
by D' (Figure S2). The value between rs2228226 (tSNP4) and 
rs2228224 (tSNPs) in UC is 0.7, providing an explanation for the 
evidence of association at rs2228226 in the absence of association at 
neighbouring SNPs. 

Found at doi:10.1371/journal.pmed.0050239.sg003 (325 KB TIF). 



ray and 



Figure S4. Quantitative Analysis of mRNA Levels b 
RT-PCR of SHH in UC Compared with Noiiiiillam, 
For the microarray data on the left, disease spccmicns are taken from 
the ascendmg colon (AC), descendmg colon (DC), and sigmoid colon 
(SC). To take into account anatomical oradienls ol HH signalling in 
Figure 3, RT-PCR (right) is perfoniu-d on samples liom SC only, and 
UC samples are siihratcgoi iscd into iioiiiiinanied (N-U and inilanu cl 

modest increase in SHH is present only with inflammation. Ihis 
result is of interest given the fact that inflammatory cells have been 
shown to express SHH. /(-Values represent Kruskal-Wallis analysis of 
HC versus UC for each anatomical location. 
Found at doi:10.137l/journal.pmed.0050239.sg004 (1.9 MB TIF). 

Figure S5. ( )iiaiititative Analysis of mRNA Levels Ijy Micro in u ind 
RT-PCR of SHH. (;LH. PTCH, and HHIP in CI) ( ipaivd with 

Noninflamed HCs 

Disease specimens are taken from the terminal ileum (TI), ascending 
colon (AC), and sigmoid colon (SC), and are subcategorised into 
noninflamed (N-I) and inflamed (I) tissues. /^-Values represent 
Kruskal-Wallis analysis of HC, CD (N-I), and CD (I) for each 

anatomical location. 



Found at doi:10.137l/journal.pmed.0050239.sg005 (1 



MB Jl 



Figure S6. Analysis of SHH, PTCH, and HHIP mRi\A Levels by 
Microarray in Colonoscopic Biopsies from the Sigmoid Colon (SC) of 
Non-IBD patients with Inflammation (I) Compared with Noninflamed 
HCs (NI) 

Found at doi:10.137l/journal.pmed.0050239.sg006 (1.7 MB JPG). 
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F^re S8. The Specificity of the N-19 HH Antibody Is Demonstrated 
with Application of the Primary Antibody, Primary Antibody + SHH 
Blocking Peptide Pre-absorbed Overnight, and with Primary Anti- 
hodv Omitted 

This figure demonstrates that the N-19 HH antibody reacts with the 
peptide but does not preclude it from cross-reacting with other 
proteins. Indeed, we would agree with other commentators [54] in 
urging caution in the use of this and other presently available HH 
antibodies for immunohistochemistry. 

Found at doi:10.137l/journal.pmed.0050239.sg008 (289 KB TIF). 

Table SI. GLIl Haplotype Frequencies in IBD, UC, CD, and HC in (A) 
Scotland, (B) Cambridge, England, and (C) Sweden 
tSNPi, rs38l7474; tSNPs, rs2228225; tSNP3, rs2228224; tSNP4, 
rs2228226. Differences between IBD, CD, UC, and HC frequencies 
are shown, calculated by test, with two-sided /(-value, OR, and 95% 
CIs. Estimated haplotype frequencies were calculated using Haplo- 
view, version 3.2. Log-likelihood /(-values (calculated on the PM/EH 
platform with 10,000 permutations and 15 degrees of freedom) are 
given for each analysis. 

Found at doi:10.137l/journal.pmed.0050239.st001 (61 KB DOC). 

Table S2. Details of the 11 GLIl SNPs Genotyped in the Scottish 
Population, Including the Four tSNPs Used for Haplotype Analysis 
X" analysis for allelic frequency and homozygosity of eai li tSNP for 
UC and CD versus HCs is presented with corres|)oii(liiig/)-yalues, OR, 
and 95% CIs. *, n. = 370 for SNPs 1-4 and 7-9; D/K, unknown. 
Found at doi;10.1371/journal.pmed.0050239.st002 (90 KB DOC). 
Text SI. Supplementary Methods 

Found at doi:10.137l/journal.pmed.0050239.sd001 (43 KB DOC). 
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Editors' Summary 

Badtground, Inflammatory bowel diseases (IBDs) are commoft disorders 
in which parts of the digestive tract become repeatedly or eontmuousiy 
inflamed. The immune system normally protects the body from entttiestt 
identifies as foreign, but in ISO it mistakenly recognizes gut tissue, and 
immune system cells accumulate in the lining of the faowet, which causes 
(nfiammation. There are two main types of tBD— Crohn's disease (CD), 
w*[fch mainly affects the small bowel, and ulcerative coiitis (UC), which 
affects only the large bowel (coion). Both types tend to run in families 
and usually develop between the a^es of 33 and 35 years. Syn^frtoms — 
induding diarrhea, abdominal cramps, and uftexplained weigtit toss — 
can be mild or severe and the disease cart develop slowly or suddenly. ;| 
There is no cure for IBD except surgical removal of the affected part of 1: 
the digestive tfact. However, drugs that modulate the immune system 

|||aj|:|||^i|i|Siiliiiiii!^^^^ by the immune s^stSii ■ f hSt stmSfsta 

||ii3||Sjii|i|||§|i||^^ in reducing symptoms. 

■'■'■'■'■'■'■'"'^ 

that harbor EHsease-BwcSciated gene variants), includinci the t8D2 locus. 
The jdentHlcation <rf tfte actaa) geie within flie IBD2 locys that is altered 
in people who are SBsefeptibte to t&D raighf provide new insights into 

liiijiiiiiiiiiliiiiiWsgiTO 
iiiliiiiiiiiiiiiipiii^sii^ 

& Do and Find? The researchers used a 



expression ol sovrr.il corTiponcnt^ of tlv- hi-dfiohofj <ii(]nnlinr| parhi'.,i> 
wrts lowor in colon s.lmpl(-^ t ik-ri Irorn p.ili' riis with UC than in wiiiploi 
l.iken Iroin hodllhy indiviciu iK =iri ill ' mice (whit h (.'xpress lidlf 
Iho nornirfl dinounl ol Cilil pr 'i in '■. -.olopod sovrro intPitiiul 
inflammation more rapidly than wild-type mice when they were treated 
vwth dextran soaium sulfate OSS), a chemical that induces acute 
(sudden) colitis. Cellular analysis revealed that myeloid cells (cells that 
sense and moaify the inflammatory response) are direct targets of the 
hedgehog signaling pathway, Furthemwre, the expression of several 
rtokincs (in particular, one called IL-23) increasec 
nn the wild-type mice after DSS 




• findings suggest that the 
'■'iijSttgiS|85SiJh inflammatory 
l|iiiiiild that reduced 
I function might be one trigger for IBD. Mots sp ili jll ih Imiiiuii 
' Mi'lii sludies identify a Gill variant that is a^,>,o( i,il'\. . illi -'J ■ ii I, jsl 
in certain north European populations), the laboratory experiments 
indicate that this GUI variant encodes a protein witii reduced activity, 
and the animal studies show that a similar reduction in GUI actMty is 
sufficient to heighten intestinal inflammatory responses. Although this 
last result needs to be confirmed in animal models of chronic col itis that 
more closely resemble human I BO, these findings suggest that drugs 
that modulate hedgehog agnaling might be usefel in #ie ttestment of 
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